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Abstract 


The wake flow around a circular cylinder at Re ~ 100 performing rotatory oscillations 
has been thoroughly discussed in the literature, mostly focusing on the modifications to the 
natural Benard-von Karman vortex street that result from the forced shedding modes locked 
to the rotatory oscillation frequency. The usual experimental and theoretical frameworks at 
these Reynolds numbers are quasi-two-dimensional, since the secondary instabilities bring¬ 
ing a three-dimensional structure to the cylinder wake flow occur only at higher Reynolds 
numbers. In the present paper we show that a three-dimensional structure can appear 
below the usual three-dimensionalization threshold, when forcing with frequencies lower 
than the natural vortex shedding frequency, at high amplitudes, as a result of a previously 
unreported mechanism: a pulsed centrifugal instability of the oscillating Stokes layer at 
the wall of the cylinder. The present numerical investigation lets us in this way propose 
a physical explanation for the turbulence-like features reported in the recent experimental 


study of D’Adamo et al. (2011 


1 Introduction 


a circular cylinder performing rotational oscillations around its axis in an infinite viscous fluid 
produces an axisymmetric pulsed boundary layer, called a Stokes layer. This is a flow susceptible 
to generate centrifugal instabilities. The linear stability problem of this flow configuration has 


been studied by Hall (19751; Seminara and Hall (19761 using asymptotic methods. A threshold 


for the appearance of three-dimensional (3D) axisymmetric instability modes was determined. 
Riley and Laurence (19761 did also stability calculations not directly on the Stokes layer problem 


but considering the modulated circular Couette flow under axisymmetric disturbances, in the 
narrow-gap limit. Later, Aouidef et al. (19941; Ern (19981; Ern and Wesfreid (1999, 20021 


considered this flow as a limit case for the stability problem of the classic geometry of two 


concentric cylinders with oscillation: the Taylor-Couette configuration (see e.g. Chandrasekhar 


(19811 for a review). In both cases, the control parameter is the Taylor number, defined as 



( 1 ) 


where i?, = ujirid/v is a Reynolds number based on the rotational angular velocity of the cylinder 
ijJi- We keep the notation of the Taylor-Couette configuration, where the subscript i stands for 
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Figure 1: Different flow states for the forcing parameters {f,A) scrutinised in D’Adamo et al. 


(20111. Visnalisations are from Thiria et al. (20061. Solid lines represent the threshold from 
global to convective instability. The dotted line indicates the threshold to centrifugal instability 
of the Stokes layer of an oscillating cylinder without a crossflow given by = 165, from|Seminara| 


and Hall (19761. Symbols ★ show the threshold for 3D centrifugal instabilities observed in the 
3D DNS discnssed in Figs. |5(b) and 5(c) Dashed region stands for states with turbnlent-like 
behavionr described in D’Adamo et al. (20111. 


inner cylinder, thns being the radius of the cylinder. In addition, TZ is the local radius of 
curvature and v the kinematic viscosity. The characteristic length scale d in the Taylor-Conette 
case is the gap between the cylinders, which fixes the scale of the wavelength of the primary 
instability. For the case studied by Seminara and Hall ( 1976[) , however, the instability occnrs 
in the inner Stokes boundary layer of thickness 5$ = y~uJoJi around the oscillating cylinder. 
They have determined analytically, nnmerically and experimentally the critical valnes for T 
associated with the onset of a Taylor-Conette-type vortex flow. Vortices evenly spaced, with a 
critical length Ac in the cylinder axial direction which is proportional to 5s, are thns developed. 


On the other hand, when a uniform flow comes across a cylinder, a prototypical 2D wake flow 
takes place for moderate free-stream Reynolds nnmbers Re = DUq/v, where D is the diameter 
of the cylinder, Uq the free-stream velocity. The well-known Benard-von Karman (BvK) vortex 


street (Benard 1908 von Karman 19111 results from the destabilisation of the steady flow in 


the wake of the cylinder and produces the periodic shedding of opposite-signed vortices with a 
freqnency /o, that occnrs above the threshold Rcc ~ 47 (see e.g. Provansal et al. (19871; Jackson 


(1987l). This flow is quasi-two dimensional up to Re < 180. In a recent work (D’Adamo et al. 


20111, we stndied experimentally the problem of the forced wake performing rotary oscillations 
at Re = 100. The rotational oscillation of the cylinder is prescribed by a forcing fnnction of 
freqnency / and amplitude 9o that can be written as 6{t) = 9 q cos{27rft), which allows the 
forcing to be uneqnivocally described using two independent non-dimensional parameters as did 


by Taneda (19781: the forcing amplitude A = ugrnax/Uo, where ugmax = Dirfdo is the maximal 


azimuthal velocity of the rotational oscillation, and the ratio f / fo- We characterised the spatial 
development of the flow and its stability properties following previons stndies by [Thiria et al.| 
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(20061; Thiria and Wesfreid (2007|. A synthesis of the case stndy is presented in Fignre 
From the analysis of power density spectra of the flow we gave a detailed description of the 
forced wake, giving insight on the energy distribution, the different frequency components, and 
in particular on a continuous spectrum observed for a high amplitude of the forcing oscillation. 
Furthermore, vortex structures revealed turbulence-like features like splitting and mixing in a 
spatial cascade pattern. A question remained concerning the physical mechanism present in the 
bifurcation that triggers such behaviour of the wake. 

We speculated on a 3D centrifugal instability to be at the origin of this sequence of transitions. 
A natural first attempt to test this idea is shown in figure where the critical Taylor number 
Tc = 165 corresponding to the instability threshold of the pure rotatory oscillating cylinder case 
without crossflow studied by Seminara and Hall ( 1976[| is id entified in the frequency-amplitude 
phase space (/, A) of the forced wake of D’Adamo et al. ( 2011[ |. This crude estimate for a 
threshold is compatible with the experimental points where the turbulent-like behaviour was 
observed (low frequencies and high amplitudes of the forcing oscillation). The purpose of the 
present paper is to characterize in detail the existence of a 3D instability and its centrifugal 
nature, using analytical estimations from the 2D flow and from 3D direct numerical simulations 
(DNS). 

It is worth mentioning that centrifugal instabilities were also reported for forced flows with 
different configurations. For transverse oscillations of a cylinder in a fluid at rest, Honji (19811, 
obtained visualizations that identified 3D structures produced by centrifugal instabilities. [Hall] 
(19841 performed a stability analysis of this configuration and gave a theoretical explanation. 


Tatsuno and Bearman (19901 investigated in detail the patterns and the structure of the flows 


that result from these instabilities. Later, Elston et al. (20041 addressed DNS calculations and 
Floquet stability analysis for this problem. 

Three dimensional instabilities in wake flows have been studied theoretically and numerically 
by Blackburn et al. (20051 where it was determined that bifurcations to three-dimensionality 
can occur from a two-dimensional time-periodic base state with space-time reflection symmetry 
for the wake of symmetrical bluff bodies. More recently for the case of the two-dimensional 
stationary flow past a rotating cylinder, Pralits et al. (20131 suggest that the stationary unstable 
three-dimensional mode could be the result of a hyperbolic instability. 


Lo Jacono et al. (20101 were interested on the role of rotationally oscillations can modify the 


three-dimensional transition in the wake of a cylinder. The frequency of oscillation was matched 
to the natural vortex-shedding frequency, /“*■ = 1, for Re = 300. They reported changes on the 
three dimensional modes from Floquet stability analysis on two-dimensional periodic flow. They 
found that the rotational oscillation dramatically suppressed mode B, even for small amplitudes 
of oscillation. Mode A was also damped, but not as significantly as mode B. For what they 
considered high rotational oscillation amplitudes, in our notation A ~ 0.66 they identify a new 
three-dimensional transition mode, which they called D mode, that shares the same symmetries 
as mode A. 

Three dimensional characteristics of forced wakes have been recently studied by |Kumar| 


et al. (20131 for the case of rotational oscillations at Re = 185 near the transition, using flow 


visualization, hot-wire anemometry and PIV. Spatial distribution of lock-on regions and its 
relationship with the forcing frequencies and amplitudes were determined. They also found 
that for certain forcing parameters (/+, A), the flow can be forced to become two-dimensional. 
Studied amplitudes were up to A = tt, a value below the threshold found in |D’Adamo et al.| 
(20111 by means of spatio-temporal spectral analysis. 

To summarise, the present work sets up a new view about 3D instabilities in wake flows, 
which have often been discussed in the case of the circular cylinder for Re > 180 as secondary 
instabilities to the BvK vortex street. We organise the paper as follows: in the next section we 
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Figure 2: a) Parameters of the problem, b) The adaptive 2D mesh around the circular cylinder 
shows the octree structure, c) Problem domain. 


describe the method used for the direct numerical simulations (DNS); results are presented in 
Sectionwhere we determine the 3D stability threshold; in Section]^ we investigate the insta¬ 
bility nature, using some concepts of centrifugal instabilities and propose therefore a reduction 
of the complex problem; lastly we elaborate our conclusions in Section]^ showing analogies with 
the Taylor-Couette problem of eccentric cylinders. 


2 Problem definition for DNS 


In order to study this problem, we performed 2D and 3D direct numerical simulations with 
Gerris free software, a parallelised tree-based adaptive solver for Navier-Stokes equations |Popinet] 
(20031. The code combines an adaptive multi-grid finite volume method and the methods of 
immersed boundary and volume of fluid (VOF). The basic equations are the incompressible 
continuity equation and Navier-Stokes equations, which can be written in terms of the velocity 
u = (u, V, w) and pressure p fields as: 


V • tt = 0 

^ + {u-\/)u = -^+iy\/^u ( 2 ) 

ot p 

The domain is spatially discretised using cubic finite volumes organised hierarchically as an 
octree. Along with the forcing problem parameters, a 2D example of the spatial discretisation 
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is given in Figure [2j The flow domain, shown in Figure 2(c) is L. 
simulations and = 20D for the spanwise direction in 3D simulations. As detailed in |Popinet] 


X Ly = 20D X lOD for 2D 


(20031 the mesh can be refined near the solid boundary and it can use vorticity gradients as an 


adaptive criterion. A cell is refined whenever: 


|V X u\Ax 
max ImI 




(3) 


where Ax is the size of the cell and ^ is a user-defined threshold which can be interpreted as 
the maximum angular deviation (caused by the local vorticity) of a particle travelling at speed 
max |m|. This adaptive criterion is represented in Figure [2(b)[ and 2(c) [ where different box sizes 
are noticeable. In order to reveal BvK vortices as well as centrifugal structures, we choose a 
minimum grid size of D/51.2 for the solid boundary and D/12.8 to define vortex regions. The ^ 
threshold is set to 0.05 for 3D simulations and to 0.01 for 2D simulations. The flow parameters 


of the simulations are defined in order to match the experimental case of D’Adamo et al. (20111: 
Cross flow velocity Uao = 1, kinematic viscosity i/ = 10“^ and cylinder diameter D = 0.1, giving 
a Reynolds number Re = 100. 


The boundary conditions are: u = 1 for x = —5D; u = 1 for y = ±5; the outflow condition 
is dv/dx = 0 and p = 0 for x = 15D; for 3D simulations, a symmetry condition is used 
for the flow at z = 20D; and at the cylinder surface, u = UgoHd where UsoUd depends on 
the forcing. As depicted in Figure rotatory oscillations are characterised by an angular 
coordinate 9{t) = 9q cos(a), where the forcing phase is a = 2'Kfft, and tangential displacements 
A = ue/{2TTff). Given /q the natural frequency of vortex shedding, the forcing frequency ff is 
written in dimensionless form as = ////o- A non-dimensional number for the amplitude of 
oscillations is obtained by comparing the maximum tangential velocity ugmax and the free flow 
velocity, A = Ugmax/Uoa- 


3 Results of the numerical simulation 


We first performed 3D DNS numerical simulations. Figure [^ shows a case with the forcing 
parameters (/“*■ = 0.75, A = 4.00). The isosurface of vorticity modulus in Figure 3(a) shows on 
one side the classic BvK wake structure synchronised with the forcing frequency. Additionally, 
a previously not reported effect is also clear: the modulation of the vorticity held along the 
direction of the cylinder axis. The two effects are depicted in Figure [3(b)[ revealing the 3D 
vortex structure around the cylinder and a well-defined wavelength A^. Moreover, Figure [^ 
shows the spatial distribution of uix along with ojz for /+ = 0.75, A = 4.00, which allows us to 
consider the symmetry properties of the observed mode. The spatio-temporal symmetry, H , of 
the two-dimensional flow is defined as: 


Huj{x,t) = KyUj{x,t + T/2) = {-uJx,ojy,-u!z){x,-y,t + T/2) 


(4) 


where Ky is a spatial reflection. For an i7-symmetric flow, from[^ the cc-vorticity changes sign 
with t —>■ t + T/2 and y —>■ —y at any fixed {x, z). This is the case for mode A, whereas for mode 
B, the sign of x-vorticity does not change. 


As studied by Blackburn et al. (20051, there are exactly three codimension-one bifurcations 
from a two-dimensional time-periodic base state to three-dimensional flow that are observable 
with variations in a single parameter. In this regard, Lo Jacono et al. (20101 showed that 


oscillatory forcing at Re = 300 leads to the appearance of a different mode (mode D) which 
has the same symmetries of mode A. Considering the symmetries observed in the present case 
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Figure 3: a) Isosurface for total vorticity modulus ||a;|j for /+ = 0.75, A = 4.00. b) Vorticity 
modulus contours for two cutting planes revealing the main flow structnres (f'^ = 0.75, A = 
4.00). 


(Fignre|^, the identified structnres are not i7-invariant, and they share the same symmetry as 
mode B. 

In what follows we thoroughly scrutinise the onset of this 3D pulsed instability. Figure |5(a)| 
shows instantaneous contours of the span-wise velocity w for a plane at y = D/2, revealing 
Taylor-Conette like vortices, with a wavelength that does not change with respect to the 
forcing amplitnde within the range 3 < A < 4. We can describe the flow with a Taylor number 
based on equation 0 considering Az as a characteristic length scale and = AUaa- We com- 
pnte the intensity of the velocity fluctuations as {w — w)'^dt/Tf, where w is the time 

average of w and Tf the forcing period Tf — I//'*'. The resnlt allows ns to identify a maximnm 
value that characterises the intensity of the 3D structure for the forcing case considered. 

An additional characterisation is possible by stndying the amplitude of this fluctuations as a 
function of the forcing parameters (/+, A). We use the 3D DNS to study the flow modifications 
for two fixed forcing frequencies /“*■ = 0.75 and /+ = 1.00 . A nseful criterion to quantify the 
intensity of 3D structnres is to follow the evolution of w. Given that the 3D structures are 
present for A = 4.00 (case depicted on Figure]^, we decrease the forcing amplitnde from this 
value until they vanish. In Fignre [5(b)[ the maxima of w^ms^ 'f^maxi are plotted against the 
Taylor nnmbers resnlting from equation Q, where the characteristic length scale Aj is found to 
be 1.16D, and the corresponding forcing amplitndes. We can appreciate, looking at the sqnare 
of the forcing amplitndes, that the 3D structnres become damped linearly as we approach a 
threshold at T = 202 for /+ = 0.75 and T = 147 for /+ = 1. The behavionr is common to 
supercritical bifurcations. 

Another scenario shows up when we follow the evolution of the intensity of 3D structnres for a 
fixed forcing amplitude. There is a range of frequencies for which the instability develops. This 
is shown in Figure 5(c) where is observed for A = 4.00 and the forcing frequency varying 
in a range 0.50 < /+ < 1.20. We can appreciate that for /+ —>■ 1.20, decreases linearly. 
On the other hand, for lower freqnencies, we observe that the 3D instability appears, with a 
finite value, for Wmax) at /+ > 0.55. 

We perform simulations for different forcing frequencies at a fixed forcing amplitude A = 4.00 
in order to characterise the evolntion of the wavelength A^. We observe in Fignre [6(a)|that A^ 
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Figure 4: For forcing parameters /+ = 0.75, A = 4.00. a) Isosurfaces for streamwise vorticity 
cUx- Wx = 3 (red), uJx = ~3 (blue). The white isosurface represents spanwise vorticity lOz = 15 
b) Top view for isosurfaces for streamwise and spanwise vorticity. 


depends on /+ following a law cx (/^)“^^^. If we assume that the “gap” size d is proportional 
to Xz, from the Taylor number definition in equation Q, where T depends on then we 
expect that high forcing frequencies produce decreasing Taylor numbers. This could explain the 
damping of 3D fluctuations for higher frequencies in Figure |5(c)[ In addition, we observe that 
the wavelength is practically invariant with respect to the amplitude for a given frequency. 
In studie s of pulsed centr i fugal instabilities. Riley and Laureii^ (1976); Carmi and Tustani- 


wskyj (19811; 


Aouidef et al. 


(19941 classified flow regimes based on a parameter 7 = 
which is the ratio of a centrifugal region length dc to the Stokes layer thickness. In our exper¬ 
iment, 7 is limited to a range between 2 and 5, it does not depend on the forcing frequency 
and dc ^ X behaves with respect to /+ as described in Figure 6 (a) where A decreases almost 
linearly as 

Even though the threshold for centrifugal instabilities determined in the Taylor-Couette 
pulsed flow is not directly applicable for a configuration with crossflow, the transformation of 
Taylor numbers based on the characteristic length dc allows an approach for our results. This 


case presents similarity with the eccentric Taylor-Couette instability problem (see e.g. Leclercq 


et al. (20131; Shu et al. (2004|; Siong (20061 and references therein). Indeed, in those problems. 


the axial wavelength of the critical perturbations is always of the same order of magnitude of 
the gap. 

Figure 6 (b) summarises the stability curves ( 7 , T) for centrifugal pulsed flow determined by 
Aouidef et al. (1994|; Seminara and (|1976[| together with the values issued from our 3D 


simulations. Two analytical curves show the solution corresponding to low values of 7 , = 

193 . 237 “^ and high values of 7 , Tc = 15 . 287 ^/^. The curves are supported with experimental 
data from Aouidef et al. (19941. On the other hand, within these reference threshold frame, 
we plotted from our results T against 7 for a fixed forcing amplitude A = 4.00, and for fixed 
forcing frequencies /+ = 0.75 and /+ = 1.00 (the same data used to construct Figure]^. We 
observe that the points are contained in the unstable region defined by the analytical curves. 
For A = 4.00, the instability develops for 0.55 < f'^ < 1.16. When /+ = 1.16, the critical point 
(7 = 3.06, T = 152) is in very good agreement with the experimental results from pure pulsed 
flows. For decreasing frequencies, T increases almost linearly regarding the estimated 7 until for 
/■*■ = 0.55 the flow stabilises with respect to centrifugal disturbances (7 = 4.4, T = 433). For a 
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(a) 



(b) 



(c) 


Figure 5: a) Contours of the transverse component of the velocity w at a cut plane y = D/2. 
Solid and dashed lines correspond to positive and negative contour values respectively, b) 
Maxima of Wrms for a fixed forcing frequency /■*" = 0.75 as a function of T resulting from 
varying forcing amplitudes, c) Maxima of Wrms for a constant forcing amplitude A = 4.00 as a 
function of the forcing frequency /+. 
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(b) 


Figure 6 : a) For a fixed forcing amplitude A = 4.00, corresponding wavelength measured in 
3D simulations, b) Taylor number versus 7 estimated from 3D simulations. Plotted values 
correspond to a fixed forcing amplitude A = 4.00 and to fixed forcing frequencies /“*■ = 0.75 and 
/+ = 1.00. Threshold points are marked in thicker symbols. The arrows represent the evolution 
to threshold from Fig. |5(b)| and Fig. |5(c)[ Our results are compared with experimental and 
theoretical data that define instability threshold in pure pulsed flows from Aouidef et al. (19941 
and Seminara and Hah] (|1976[| 


fixed frequency /“*■ = 0.75, 7 = 3.61, the flow destabilises at T = 202 and with increasing forcing 
amplitudes, T eventually reaches the previous set of points at A = 4.00. The same behaviour is 
found for the fixed forcing frequency /“*■ = 1.00, where the flow is unstable from T = 147. 

We suggest that the centrifugal instability that develops in the forced wake can be thus 
considered in the context of pure rotatory pulsed oscillations. Nevertheless, the natural BvK 
dynamics plays an important role as the first bifurcation depends on the distance of the forcing 
state space parameters (/, A, Figure[^ to the resonance centred at (/+ = 1.00, A = 0.00). This 
fact could explain that at < 0.55 the centrifugal instability is not strong enough even when 
T is high. Conversely, for (/+ = 1.16, A = 4.00) the length dc is significantly smaller and T 
decreases to the values predicted by the pure pulsed flow threshold. 

We bring a quantitative picture of these ideas in the remainder of the paper, starting with 
a brief review of the criterion for centrifugal instability. 


4 Centrifugal instability 


The necessary condition for a 3D centrifugal instability in flows with curved streamlines is given 


by Rayleigh (19161 criterion for inviscid flow (see e.g Drazin and Reid 19811 which can be 
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Figure 7: In non-forced wake, instantaneous streamlines and (j) negative contours reveal poten¬ 
tially centrifugal instabilities regions. Regarding curvature and velocity modulus, two distinct 
regions associated with Gortler and Taylor-Couette mechanisms are marked. However, 
—(j) is not strong enough to overcome the stabilising effect of the viscosity and to produce 3D 
instabilities for Re = 100. 


written for flows such as the Taylor-Couette flow in terms of the Rayleigh discriminant 


, 2V fV dV\ 


(5) 


where V{r) is the 2D velocity of an orthoradial base flow field. 3D perturbations to this flow 
field are amplified if < 0, which translates the fact that the perturbed pressure held does 
not balance the centrifugal force, leading to flow instability. For a general profile V{r), the 
flow field can be subdivided in regions of different stability depending on the sign of it 

will be unstable in the region where ^(r) < 0 and stable when (f){r) > 0. More generally, for 
other geometries described by a vorticity held uiz, the Rayleigh discriminant can be written as 
(j) = V{r)uiz{r)/TZ ( Leblanc and Cambon[ 1998, Sipp and Jacquiii] 20001, where the local radius 
of curvature of the streamlines TZ is defined by 


7 ^ = 




{Atp){u\7u) 


( 6 ) 


Beaudoin et al. (2004| made use of these expressions in order to identify potential instability 


regions in a backward-facing step flow and characterise the 3D global instability. In the present 
case, we will see that the study of the local Rayleigh discriminant is a useful tool to predict the 
centrifugal stability of the forced cylinder wake problem, which lacks of symmetry simplifications. 
In Figure]^ the instantaneous flow streamlines along with the Rayleigh discriminant 4>{x,y) are 
represented for a non-forced flow around a cylinder at Re = 100, where the flow produces the 
Benard-von Karman vortex shedding. Two distinct regions of potential centrifugal instability 
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Figure 8: Rayleigh discriminant negative values contours, (^ < 0 for /+ = 0.75 and A = 4.00, for 
three different times or oscillation phases within half of the forcing period Tj. The left picture 
shows (j) spatial distribution when a vortex is shed below x = 0 axis, when the forcing phase 
a ~ 1 . Middle (a ~ 2) and right (a ~ 3) pictures allow to appreciate the temporal and spatial 
evolution of (j). Symbol ir stands for the position Xmax of |^|max- 


exist: one near the stagnation point, where a concave streamline constitutes a Gortler-like 
geometry (Saric, 19941; and another one in the near wake side, where the curvature of the 


streamlines around the vortex formation region corresponds to a Taylor-Couette geometry. At 
Re = 100, nevertheless, viscosity prevents the development of 3D instabilities, which never 
appear for the case shown in Figure When the rotational oscillatory forcing is applied, 
negative values of (p appear mostly in regions close to the cylinder. In what follows we define 
the characteristic length scale dc of regions potentially unstable giving a local Rayleigh criterion 
to analyse the stability properties of the forced wak^ 

Despite the flow complexity, it is possible to reduce the problem to investigate solely the 
centrifugal instability of the 2D base flow and its relationship with the forcing parameters. We 
calculate the Rayleigh discriminant (j){x, t) for each forcing parameter from the streamlines of 
the flow at Re = 100. In Figure]^ we present three snapshots of (/>(x, a) at /+ = 0.75, A = 4.00 
and Re = 100, where a € [0, 27r] is the forcing phase. For other forcing parameters we obtain the 
same qualitative features than what we describe for Figure We observe that a “corona”-like 
region appears around the cylinder with negative values of p. Figure]^ a) shows the phase when 
(p is the most negative, where we can expect the strongest possible centrifugal instability with 
the highest growth rate ( |Bayly (1988)). The location Xmaxj where the instability can be the 
most strong locally, is given by |(/)(xmax)| = |<)>max|- 

Figure |^b) and [^c) describe the evolution of p{x, a) for forcing phases that correspond to 
the mean and the minimum values, where the flow is less receptive to the instability. We observe 
that p{x, t) is x-symmetric regarding the forcing phase, <p{x, y, a) = <p{x, —y, a + tt). In Figure 
we show the variation of the lift coefficient cl = L/^pU^, being L the resulting lift force, 

^It should be noted that viscous effects have to be considered to determine the actual stability criterion. This 
could be done using dc in the definition of the Taylor number T of equation^ 


11 





















u 


Figure 9: |0|maxj red circle symbols, is the maximum of the Rayleigh discriminant modulus in 
function of the forcing phase a, |(/)|max = |'i^’(a;max,a)!- The lift coefficient cl, green triangle 
symbols, gives a reference for the forcing phase a. The local radius of curvature value TZ, brown 
diamond symbols, at Smax is close to the cylinder radius, its sign changes with the shedding 
cycle. Blue square symbols ■ mark three forcing phases depicted in Fig. 


which is correlated with the phase reference a. The local radius TZ calculated from Eq. ([^ in 
^max is also represented, we can see that its modulus is close to the value of the cylinder radius 
as the curved streamlines of the forced flow approach the cylinder. The sign of the local radius 
accompanies the changes due to the oscillation. We present in the same Figure a curve for the 
evolution of the maximum of the Rayleigh discriminant modulus t)|max during a forcing 
period. We indicate over this curve the three values of |0|max that lead to the construction of 
Figure [8| 

As we have already pointed out, we can extract a convenient length scale in order to adapt 
our problem to the pure centrifugal instability framework, allowing us to compare our results 
with previous works. If we consider the forcing phase that corresponds to |(^|max, around a = 4 
in Figure the 2D flow streamlines are depicted by Figure [T^a). Using equation the 
Rayleigh discriminant is obtained and presented in Figure [T^b). It is worth mentioning that 
a y-symmetric field is retrieved for a ~ 1 that corresponds to the other maximum of \(j)\. 
Image processing is used in order to extract a length scale from a contour plot of the Rayleigh 
discriminant obtained from equation a s shown in Figure [T^c) (see Appendix for details). 
The mean radius represented in Figure]lO|(c) determines the length scale dc related to the size 
of the unstable region for 2D flow. It is shown in Figure ll(a)| for different forcing frequencies 
at a fixed forcing amplitude A = 4.00 together with the size A^/2 of the centrifugal rolls that 
develops in the 3D flow. Both dc and \z/2 follow the same trend, supporting the idea 

of the pulsed Stokes layer. Their ratio, around a value of 3, is plotted in Figure |ll(b)[ Given 
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c) 

Figure 10: For (/+ = 0.75, A = 4.00) at the forcing phase a^ax corresponding to ^max = 
(^(x, amax): a) Streamlines, b) Regions of negative Rayleigh discriminant —cj). The spatial 
minimum of (j) is placed at a position identified with the symbol tV. c) </> regions after image 
processing. The corresponding mean radius for (/) is represented and the resulting measure dc is 
indicated. 
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Figure 11: For a fixed forcing amplitude A = 4.00, variation of the estimated centrifugal insta¬ 
bility measure dc at a ~ 4 from 2D DNS and variation of the wavelength associated with the 
size of the rolls in Figure 5(a) b) Ratio between the characteristic length dc for the centrifugal 
instability determined from 2D simulations and Az/2 issued from 3D simulations. Th thick black 
line represents the mean value of Az/2. 


that the flow is under non-stationary forcing, the rolls are formed periodically symmetric with 
respect to the x-axis. Besides, dc has been determined for a particular phase a ~ 4, where 
the instability is most intense, but the centrifugal instability region changes its size. These 
arguments may explain the difficulty for estimating dc and the the scale difference between Az/2 
and dc- Nevertheless, we can observe that the main behaviour is shared between Az and dc, 
therefore, a centrifugal instability region observed in 2D simulations is in agreement with the 
3D instability that develops in 3D DNS. 

5 Conclusions 

The present work gives a a new view about 3D instabilities in wake flows. In the context of forced 
wakes at moderate Reynolds numbers, we found a new transition that leads to the formation 
of three dimensional structures. The instability shares aspects that were previously studied 
for centrifugal pulsed flows. Taylor Couette-like vortices develop from a definite threshold of 
forcing parameters (/, ^4) and these structures are modified by the incoming flow. For this 
complex instability, 2D evaluation of the Rayleigh discriminant may give a fast criterion to 
determine whether a wake flow becomes three dimensional or not. We found from streamline 
shapes and the spatial distribution of (j) that the problem shares some analogy in relation to 
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eccentric Taylor Couette flows. 

As 2D forcing in wakes may indeed trigger 3D structures, this behaviour must be taken into 
account in flow control schemes. Streamlines which result too much “bent“ by forcing in wakes 
can make evident strong negative values of the Rayleigh discriminant cj) and thus the possibility 
of a centrifugal instability. 

On the other hand, this simple problem can offer an interesting benchmark to study instabilities 
and transition to turbulence from oscillatory rotation. 


6 Appendix: Determination of the centrifugal instability 
region length. 

The choice of a characteristic length of the centrifugal instability region from the Rayleigh 
discriminant scalar fields is not straightforward as we observe Fig. we choose to select the 
forcing phase that corresponds to the minimum value of (j), the most unstable state. Figure [lOi b) 
presents such state but the (p scalar held needs to be more clear in order to extract a length dc- 
Simple image processing functions, erosion and dilation, are applied successively to the scalar 
held in order to obtain Fig [Tolc), where a clear shape is noticed. We found that such shape 
has an aspect that resembles an eccentric cylinder gap. Therefore, we choose as a characteristic 
length the mean radius of this gap dc = l/(27r) rdip — D/2, with r the shape radius varying 
with the angular coordinate cp. 


6.1 Convergence analysis for DNS 


In order to ensure that the results do not depend on the size of the domain we chose, we 
performed a convergence analysis for the 2D case. Given that the domain size of the reference 
study is La; X Lj, = 20L> x lOD, we label it as LI.00. As we selected larger domains which scale as 
[1.25; 1.50; 1.75; 2.00] the reference study length, we label them LI.25, LI.50, LI.75 and L2.00. 
For these scaling lengths, we plotted mean flow profiles for the streamwise component of the 
velocity Um at three different x positions x = 0.5D, D, 2D in Fig. 12(a) Figure [T2(b)| presents 
for the same direction, fluctuations intensity Urms profiles for the same different x positions. 
We observe that changing the domain size does not modify the flow dynamics. Lift coefficient 
is also calculated for each case as it is presented in Fig{^ We also observe good agreement 
between the different scaling domains. 


Data Accessibility 


In order to reproduce all the calculations included in this paper, Gerris(Popinet 20031 is avail¬ 
able free of charge under the Free Software GPL license and our code files can be downloaded 
from this URL, 
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(a) 



(b) 

Figure 12: For the forcing case f~^ — 1.00 and A = 3.50 we test the reference scale LI.00 versus 
larger scales up to L2.00 . (a) Mean flow profiles Um at x = 0.5D, x = D and x = 2D for the 
different scales tested, (b) Fluctuations intensity Urms for flow profiles at the same locations. 



Figure 13: For the forcing case = 1.00 and A = 3.50, the lift coefficient Cl evolution during 
two forcing periods corresponding to five different scale domains. 
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